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SUMMARY 


The compressive Buckling stress of outstanding flanges reinforced 
By BulBs was computed By the strain-energy method for flanges having 
10 shapes and a range of lengths. The results were checked for some 
cases By computations Based on a differential-equation method. The edge 
of the flange opposite the bulB was considered claused, and the loaded 
ends were considered simply supported. The results were analyzed to 
determine which shape of flange gave the greatest support to the struc- 
ture to which it was attached. It was found that the flange capable of 
giving the most support without Buckling had on over-all flange width 
of 3.4 ^Aj^ (where Ap is the cross-sectional area of the flange). 

This flange was somewhat slimmer than the one previously found Best for 
the case of a simply supported flange root. 


INTRODUCTION 


Flanges reinforced by BulBs are widely used to reinforce the 
stressed skin in semimonocoque aircraft structures . They perform a two- 
fold function in the structure By (1) maintaining the structural shape 
and (2) carrying a substantial portion of the externally applied load. 

The outstanding flanges in such structures have a tendency to fail 
under compressive load by a lateral displacement of the outstanding 
portions of the flange with respect to the rest of the structure. As 
is pointed out in reference 1, such a failure is intermediate Between 
torsional instability in which no cross-sectional distortion of the 
flange occurs (considered By Goodier in reference 2, Kappus in refer- 
ence 3, and Lundquist and Fligg in reference 4) and local instability 
of component elements of the structure in which the lines joining 
component plates remain straight (considered By Timoshenko in refer- 
ence 5 > Kroll in reference 6 , and Lundquist, Stowell, and Schuette in 
reference 7)- 
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The lateral instability of outstanding flanges under compressive 
load, for the case where the flange has a reinforcing bulb at its outer 
edge, has received considerable study. Windenburg (reference 8) consid— „ 

ered both local and primary instability of flanges having a reinforcing 
bulb symmetrically placed on both sides of the flange web. He consid- 
ered the case of simple support of the flange by the structure to which 
it is attached. Chwalla (reference 9) considered the case where the 
structure gives clamped support to the flange and presented numerical 
results for certain specific cases. 

A survey (reference 1) of aluminum-alloy extrusions showed that 
the shape of most reinforcing bulbs is either rectangular with a rounded 
end and a fillet at the junction of the bulb and flange or circular with 
a fillet at the junction of the bulb and the flange. In nearly all 
cases the bulb was entirely on one side of the flange, thus resulting 
in an unsymmetrical distribution of material at the outer edge of the 
flange. It was found that the range of cross-sectional shapes used 
could be approximately covered by 54 specific flange shapes. 

In reference 1, numerical solutions were given in dimensionless 
form for the lateral instability of these different flanges as a 
function of the flange length. The edge condition along the line of 
attachment to the structure was taken as the extreme of simple support. 

In the present paper, the case of clamped support will be treated and, 
of the 5 k cross sections analyzed in reference 1, only a representative 
group of 10 will be considered. 

This work was conducted at the National Bureau of Standards under 
the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 




cross-sectional area of bulb; bulb includes area cross- 
hatched in figure 1 


Ajp 

As 

b 

*1 


cross-sectional area of flange 

cross-sectional area of sheet in sheet— stringer structure 
flange width from base to bulb center 

width of rectangular— type bulb from free end to center line 
of web 
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flexural rigidity of wet 


Young’s modulus 


12(l - v 2 ) 


shear modulus 


2(1 + V) 


Poisson’s ratio; 0.3 


moment of inertia 


moment of inertia of bulb about center line of web 
moment of inertia of flange about its base 
torsion constant of bulb 
flange length 
flange thickness 

thickness of rectangular-type bulb 

fillet radius of rectangular— type bulb, taken equal to tj 

fillet radius of circular-type bulb 

bulb radius of circular-type bulb 

lateral displacement of flange 

compressive end stress 

critical stress for instability of flange 
critical stress for sheeb-stringer structure 

critical stress for instability of an infinitely long 
flange 

angle of twist of bulb with respect to web 
coordinate along flange length 


coordinate along flange width 
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SHAPE OF REINFORCED FLANGES 


The shapes of the reinforced flanges are shown in figure 1 and 
their properties and dimensions are given in table 1. Five of the 
flanges had essentially rectangular reinforcing bulbs, while, for the 
remaining five, the bulbs were essentially circular. The first column 
in table 1 gives the flange number. The second column gives the number 
of the corresponding flange considered in reference 1. The third 
column gives the thickness of the flange web, and columns (4) to (6) 
give dimensionless ratios describing the shape of the cross section. 


ANALYSIS 


The flanges considered in this report are clamped at the root edge 
and get lateral and torsional support at the outer edge from the 
reinforcing bulb. For such a flange, no bending is possible in the 
plane of the web and, therefore, any bending which does occur will have 
the web center line as a neutral axis. This fact will allow a conven- 
tional treatment for the unsymmetrical reinforcing bulbs (fig. 1) 
considered in this report. The buckling stress of the flange will be 
confuted as that of a plate clamped at one edge and restrained at the 
other edge by a beam. This beam will have a bending stiffness EI^ 
equal to that of the cross-hatched area of figure 1 about the center 
line of the web, a torsional stiffness GJ-j equal to that of the whole 
flange (fig. 1) less that of the web, and an area A-^ equal to that of' 
the cross-hatched area of figure 1. Values of 1^, JVp, and Afj are 
given in columns (7)* (8), and (9) of table 1. The total area of the 
flange Ap is given in column (10) of table 1. The flexural rigidity D 
is given as D/E in column (ll) of table 1. 


Plate Theory 

An expression for the compressive buckling stress was first 
obtained by integrating the differential equation of the deflection 
surface to find the magnitude of the compressive stress necessary to 
keep the plate in a slightly buckled shape. The derivation of this 
method appears in the appendix. Since the integration method appeared 
excessively cumbersome, a strain— energy method was developed, using the 
technique given on pages 325 to 327 of reference 5 as follows: 
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The work of the external forces on the plate and team (see p. 325, 
reference 5) is 



The strain energy of the team and plate, neglecting shear in the 
team is (see pp. 3^5> 24, and 265, reference 5) 


(D/2) 






\'~ 2(1 - V) 


^2 W ^2w 


dx 2 8y 2 


^2w \ 2 

f&dy) _ 


dx dy + 



dx 


0J T 

♦-2 1 


PL 



dx = Ii 


(2) 


where 



In order to obtain the results in dimensionless form it was 
convenient to set v = 0.3, G = E/2.6 and to work with Ip’, defined 
by 


I I’ = Il/E = (D/2E) 


’t al r 


0 J 0 


'^S + ^»Y-14 


c)2w ^2 V / ^2 V \ 2 

8x 2 By 2 \(8x8y/ J 


*b 

2 


’ L (Wf fW 

\^V y =b 5 ' 2 Jo V to /y=b 


f dx: dy - 


00 


A function allowing a single buckle in the x-direction. 


w = ( sin 


f Vai/ ♦ a ^ 3 * a 3 / . . . V n+1 ) 


( 5 ) 


where a^, ag, . . ., a^ are arbitrary constants, was chosen to 
represent the deflection surface of the plate. 
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P* 


The critical stress is given "by the expression (reference 5, 

326 ) 




or, using equation (4), 


X 2 a cr ~ I 1 



( 6 ) 


(7) 


The values of the constants a-,, a^, ... in equation (5) are 

determined from the condition that they give the lowest possible 
"buckling stress <j . For this to "be true they must satisfy the set of 

Ci 

simultaneous equations (reference 5, p. 326) 



g cr _ 
E da.^ 


dll* + ^cr d*2 = 0 

da 2 E 


( 8 ) 


t • t • • 

Solutions of equations (8) different from zero may he obtained 
only if the determinant of the equations is zero. 


The expressions for It' and I 2 may be obtained by substituting 
equation (5) in equations ( 4 ) and (1), respectively. The critical 
buckling stress may then be obtained by substituting the expressions 
obtained for I 1 ‘ and I 2 in equation (8), setting the determinant 
of the resulting equations equal to zero, and solving for the lowest 
values of a cr /®* 


Convergence 

To determine the minimum number of arbitrary constants necessary 
in equation (5) for accurate determination of the buckling stress, the 
buckling-stress ratio of flange 5 with lA/A^ = 5-396 was computed 
three times, using one, two, and three arbitrary constants in equa-- 
tion (5), successively. The buckling-stress values according to the 
strain— energy method were as follows : 
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Number of arbitrary constants 

Values of ff cr /E obtained 

in deflection equation 

1 

0.0308 

2 

.0282 

3 

.0233 


The convergence seemed slow. To obtain additional information on this 
point, the buckling stress was also confuted by the integration method 
with the result •^cr/E ~ 0.0233* 

It was concluded that in this case the equation 

w = (sin + a 2 y3 + a^) ( 9 ) 

has the minimum number of arbitrary constants necessary to adequately 
represent the deflection surface. 

Buckling stress— ratio values for one or more length— ratio values 
were computed for each of the 10 flange shapes by the integration 
method given in the appendix for comparison with the corresponding 
buckling stress obtained by the strain— energy method using equation (9). 
Stress— ratio values obtained by both methods are given in table 2. In 
all cases the results of the two methods agree within the accuracy of 
the numerical computations. It was concluded that equation (9) would 
be sufficiently accurate and the remaining computations were done by 
the strain— energy method using equation (9). 


Contribution of Bulb to Critical Stress 

The contributions of the bending stiffness and torsional stiffness, 
respectively, to the buckling stress of flange 4 at = 18.91 were 

determined by recomputing the buckling stress ratio twice, first 
omitting the bending— stiffness term, and then omitting the torsional- 
rigidity term. It was found that omission of the bending— stiffness term 
reduced the buckling stress— ratio value from 0.0147 to 0 . 0086 , or 
41 percent, and that omission of the torsional— stiffness term reduced 
it from 0.0147 to 0.0134, or only 9 percent. 


Multiple Roots 

In the strain— energy method, the critical stress is found from the 
condition that the determinant of the coefficients of three simultaneous 
equations shall be zero. This determinant, when expanded, is a cubic 
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equation, and as a result more than one solution is possible. In order 
to make sure that the solution obtained was that corresponding to the 
lowest buckling load, a plot was made, in some cases, of the value of 
the determinant as a function of the assumed critical stress, starting 
from zero stress. The first intersection with the zero axis was taken 
as the lowest critical stress. This method was used only in those 
cases where there was some doubt as to whether the stress obtained 
actually corresponded to the lowest critical stress. In the other 
cases, values of the determinant were Qbtained only in the neighborhood 
of the estimated zero point. 


RESULTS 

Critical Stress 


The buckling stress ratio cr cr /E is plotted in figure 2 against 
the length ratio L/\jAp for the 10 flange shapes considered. The 
ratio L //% was cnosen as the abscissa since it has a fixed value for 
all geometrically similar specimens and it also has a fixed value for 
all geometrically different specimens of a given length having the same 
cross-sectional, area. Comparison of ordinates for a given abscissa 
will therefore show the effect of changing the distribution of material 
in the flange. 

The buckling stress ratios in figure 2 decrease rapidly with 
increasing flange length, go through a minimum at a flange length lower 
than that usually encountered in aircraft structures, and rise rapidly 
for greater flange lengths. Stubby flanges are more stable than wide, 
thin flanges, and, except in the case of the two narrowest flanges, 
flanges reinforced by rectangular bulbs are more stable than those 
reinforced by round bulbs. 

Flanges as long as those commonly used in aircraft structures will 
buckle longitudinally into a number of half sine waves. In the case of 
flanges with simply supported loaded ends the number of half waves into 
which the flange buckles and the buckling stress for any length may be 
obtained by plotting a family of curves, using the ordinates of the 
original stress— ratio, length— ratio curve, and multiples 1, 2, 3, and 
so forth of the corresponding abscissas, successively, as in figure 3- 
The transition from a buckle pattern of one half wave to that of two 
half waves occurs at the intersection of the first and second curves so 
drawn, from two to three half waves at the intersection of the second 
and third curves, and so forth. The critical stress ratio for any 
length ratio may be obtained from the curve having the lowest ordinate 
at that length ratio. With increasing length ratio the buckling stress 
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ratio approaches a constant minimum value a cr m ^ n / E * Deviations from 

this minimum value are generally small in the range of flange 
lengths L > 20^A^T. 


Transverse Deflections of Weh Center Line 

An expression for the relative transverse deflection of the weh 
center line of each of six flanges subjected to buckling stress was 
obtained as follows: Two of equations ( 8 ) were arbitrarily selected 

and the buckling stress ratio was substituted for <J cr /E. The coeffi- 
cient &3 was then factored out and the resulting equations were 
solved simultaneously for ap/a^ and /a^ . Using equation ( 9 ), the 

relative values of the deflections were confuted. 

These deflections are plotted in figure 4, for three flanges with 
rectangular bulbs, and in figure 5 > for three flanges with circular 
bulbs. In figures 4 and 5 , the ordinate is the ratio of the deflection 
at any point on the web center line to the maximum deflection, and the 
abscissa is the ratio of the distance of the point from the flange base 
to the total flange width at the web center line. The flanges are 
drawn to scales of equal area on the figure. The flange— deflection 
curve starts at zero slope from the built-in edge, curves sharply, and 
straightens at the bulb. Flanges having relatively large bulbs tend to 
have s— shaped deflection curves of the web center line. 

The effect of flange length on the shape of the transverse deflec- 
tion curve is shown in figure 6 , using flahge 4 as an example. A range 
of flange lengths from L /\Jk$ = 9*^-5 bo 27. 01 was considered, since it 
covers the range of lengths near the minimum buckling stress value as 
shown by figure 3* The web center line tends to be s— shaped for low 
values of L/^Ap, indicating that the effective stiffness of the bulb 
increases rapidly with decreasing flange length. 


EULER INSTABILITY OF SHEET-STRINGER STRUCTURES 


A method for determining the effectiveness of the flanges in 
preventing instability by Euler column buckling parallel to the plane 
of the flange was outlined in reference 1. It was concluded that, if 
the structure was not to have Euler column instability before a 
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stress 0 crg , the largest area of sheet Ag which a flange of moment 
of inertia Ip about its base would stabilize was given by 



( 10 ) 


The most efficient flange for a given stress ratio cs crs j E and 
length ratio l/^Ap will be that having the highest value of Ag/Ap. 

Equation (10) shows that this corresponds to having the largest value 
of Ip/Ap > that is, radius of gyration relative to the base of the 
flange for a given cross-sectional area. This conclusion assumes that 
no other instability occurs before column failure. It rules out flanges 
with large values of Ip/Ap 2 which fail by any other instability at 

relatively low stress. Inspection of column (l4) of table 1 shows that 
flange 4 has the highest value of Ip/Ajp 2 and is closely followed by 
flanges 9, 6, 5, and 10. Inspection of figure 3 shows that, of this 
group of five flanges, flange 4 also has the highest local buckling 
strength c cr /E for all values of l/\/Ap and that this strength is 
well into the plastic range since it is never less than 0.0147- 
Flange 4, therefore, is the best of the 10 flanges considered in this 
report . 

The area of sheet Ag, which flanges having the shape of flange 4 
can support, is shown in figure 7 for a range of values of a crs (the 
Euler column instability stress of the structure) from <Xcr s /E = 0.004 
to 0.015. The ordinate in this figure is Ag/Ap while the abscissa 

is L/y'Aj?. The area of sheet Ag which can be stabilized decreases 
sharply as either ff C r s /E °r Iy^A^ increases. At high stress, 

a CT / E = 0.012, and reasonable length, L/^/3Sp = 40, Aq is slightly 
less than twice Ap. 

COMPARISON OF PRESENT FLANGES WITH SIMPLY SUPPORTED FLANGES 


It is of interest to compare the performance of the present flanges 
clamped at the root with the performance given in reference 1 for 
similar flanges simply supported at the root. 
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In reference 1 it was found that the most effective flanges had 
relatively compact cross sections. The best over— all flange width 
ranged from 1.9 yfA^ to 2.6 \jAj? and hulb width from 0-7 /Ay to 1.2 \p^. 
In the present report the best flange has a somewhat slimmer appearance, 
with a flange width of 3.h and a bulb width of 1.2 ^jA^. 

National Bureau of Standards 

Washington, D. C., October 17 > 19^7 


12 


NACA TN 1985 


APPENDIX 


SOLUTION BY INTEGRATION METHOD 


A solution for the "buckling strength of a plate built in at one 
edge and restrained by an elastic beam at the opposite edge may be 
obtained as follows: 


The equation of the plate surface is (reference 5j P- 324) 
6^w + ^ S^w + c)^w _ _ q xt d^w 

cix 2 oy 2 D 8x 2 

At the edge y = b (reference 5, p. 343 ) 


(Al) 


GJ q 


83w 


_/8 2 w 8 2 w \ 

+ V/ = 8x 2 dy 


(A2) 


for torsional edge condition, and at the edge y = b (reference 5, 
p. 346) 

A ■ * r S Jij 


EI 




Vx 8x 2 


(A3) 


for lateral— force edge condition. At the edge y = 0 (reference 5, 
P- 344) 


w = 0 
8w 


8 y 


= 0 


( A4) 


Following pages 341 and 337 of reference 5 and taking m = 1, 


JtX 

w = sin — 
L 


A(cos Py — cosh ay) + B^sin py - ^ sinh ajj 


(A5) 





NACA TN 1985 


13 


where 



and A and B must satisfy the conditions of equations (A2) and (A3), 
These give two simultaneous equations for A and B which will only 
have nonzero values for a particular value of o x . 

The two conditional equations, with expressions for the partial 
derivatives obtained from equation ( A5 ) substituted, sin 252-j 

factored out, and y set equal to b, are: 

( o ? vji 2 vjt 2 \ 

Alp^ cos pb + a cosh ab + — cos pb — cosh ab ) + 

\ L 2 L 2 / 


B^p 2 sin pb + ap sinh ab + sin Pb — - sinh ab 
\ L 2 aL 2 j 


GJrr 


_.2q 

A( £ sin pb 


sinh ab^ + B (^~ cos pb 


it 2 p 


cosh ab 


(A6a) 


and 

Hf ♦ V * | 


A(cos pb — cosh ab) + B^sln pb — — sinh ab 


{s 


A^— p3 sin Pb + a ? sinh ab ) + B^p^ cos pb + a 2 p cosh ab) + 

* 2 n 2 

(2 — v)Av -o(— P sin pb — a sinh ab) + (2 — v)B-^(P cos pb — P cosh ab) 

L 2 L 2 

(A6b) 


Solving equations (A6a) and. (A6t>) for A/B gives 
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For a given flange geometry the quantities A^, b, L, J T , 
D/E, and t are known. The value of v is taken as 0.3. 

For a given value of a x / E there can be computed 


+ ^(10.92) ■£* ■ r 


L 2 V E 


t 2 L 2 


Pb =\/-rt 2 777 + \/^r * 2 (10.92) 


VE 




= s 


CSJrpJt^ ^ 1)2 

L 2 Eb 2 " 2.6b^ L 2 J 


b 2 /L 2 = c 


t3 


Eb3 10.92b 3 


= h 


~ EI $ + Vx ^)| * -” 1 * jr * x 2 ft ft - 


Let 
U = 


js(cos s — cosh r) — h|s 2 sin s + rs sinh r + 2.96lc(sin s— p- sinh r)J 
j(s sin s + r sinh r) + h £s 2 cos s + r 2 cosh r + 2. 961c (cos s— cosh r )J 


hs£s 2 cos s + r 2 cosh r + l6.78c(cos s— cosh r)J — p^sin s— p sinh r) 
h£s 3 sin s — r 3 sinh r + l6.78c(s sin s + r sinh r)^+p(cos s— cosh r) 


The intersection of U and 1 plotted as functions of crx/E* 
defines <r cr /E. 
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TABLE 2.- BUCKLING- STRESS-RATIO VALUES OBTAINED BY 
BOTH STRAIN-ENERGY METHOD AND INTEGRATION METHOD 


Flange 

l/\Z^f 

“cr/® 

Obtained, by- 
strain— energy method 

Obtained by 
integration method 


12.47 

0. 0244 

0.0244 

1 

19.96 

.0185 

.0184 

O 

5.23 

• 0570 

• 0570 


17.43 

.2207 

.2207 

3 

11.62 

.0236 

.0236 

4 

18.91 

.0147 

.0147 


5.40 

.0233 

.0233 

5 

29.98 

.0245 

.0245 


7.18 

.00648 

.00648 

6 

8.98 

•006l4 

.00614 


5.88 

. 0666 

.0666 

7 

7.35 

.0642 

.0642 

Q 

4.91 

• 0333 

• 0333 

O 

9.83 

.0187 

.0187 

9 

12.92 

. 00880 

.00880 

10 

5.76 

.0270 

.0270 

14.41 

.0113 

.0113 








subjected to buckling stress. 



subjected to buckling stress. 



Figure 6.- Amplitude of transverse deflections at web center line for flange 4 at buckling stress. 






